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Abstract 

We investigate the task of controlling ensembles of initial and terminal state vectors of parameter- 
dependent linear systems by applying parameter-independent open loop controls. Necessary, as 
well as sufficient, conditions for ensemble controllability are established, using tools from complex 
approximation theory. For real analytic families of linear systems it is shown that ensemble con¬ 
trollability holds only for systems with at most two independent parameters. We apply the results 
to networks of linear systems and address the question of open-loop robust synchronization. 
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1 Introduction 

Driven by recent engineering applications the task of controlling ensembles of system by open loop 
controls has gained considerable attraction. The motivation for this study originates, for instance, 
from quantum control [13], the control of spatially invariant systems such as platoons [U [3], the 
control of partial differential equations [3|, and flocks of systems |2]. 

From a functional analytic point of view, the problem of ensemble controllability is equivalent to 
approximate controllability of infinite-dimensional systems defined on Banach or Hilbert spaces |9]. 
Standard characterizations of approximate controllability in Hilbert spaces can be found in mm- 
However, the results in [3 [ID are, except for very special cases, not easily applicable for ensemble 
control as they depend on the existence of a Riesz basis of eigenvectors. A new function theory 
approach to uniform ensemble controllability has been developed in [S|, using classical approximation 
theoretic results, such as the Stone-Weierstrass Theorem |6| and Mergelyan’s Theorem |18| . 

In this paper, recent results on ensemble control by PESKini are extended in several directions. 
The necessary and sufficient conditions for ensemble control in |9] are extended to finite unions of 
disjoint parameter intervals. Somewhat surprisingly, we show that ensemble controllability of real 
analytic families of systems holds only if the parameter space is at most two-dimensional. New 
explicit characterizations of uniform ensemble controllability are derived for a special class of one- 
parameter families of systems. Finally, we discuss an application to robust synchronization for circular 
interconencted homogeneous networks. 
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2 Polynomial characterizations of uniform ensemble controllability 

We consider parameter-dependent linear time-invariant systems, described in continuous-time by 


9) + B{e)u{t) 

[x(O,0) = O, 


( 1 ) 


and in discrete-time, by 


S(0): 


xt+i{9) = A{e)xt{6) + B{9)ut 
xo{9) = 0. 


(2) 


Since we want explore reachability properties of these systems, the initial state is taken to be zero. 
Note that in the continuous-time case there is no difference between reachability and controllability. 
In the sequel we will only use the term controllability, with the caveat that for discrete-time systems 
this has to be interpreted as reachability. Recall that a system © , or m, is reachable for a fixed 
parameter 9 if and only if the Kalman matrix has full rank, i.e. 


rank(R(0) A{9)B{9) ■■■ A{9Y-^B{9)) = n. 

Throughout the paper we assume that the parameter 9 varies in a compact subset P C and the 
system matrices A{9) € are continuous functions in 9, i.e. A € C'(P,M”'^"') and 

B € (^(P, We refer to such parameter-dependent families of systems and initial/final states 

as an ensemble. Throughout this paper we assume that P is a nonempty compact subset of with 
int P = P. This implies that the dimension of P is well defined and satisfies dimP = d. It excludes 
the well-understood case of parallel interconnected linear systems, where P is a finite set. 


In this paper we address the following open loop control task for parametric families and ©. 
Let 

>fiT,9,u)= [ e^(^)(^-")R(0)u(s)ds 

Jo 

and 


T-l 

^{T,9,u) = J2M0)’^B{9)uT-i-k, 

k=0 

denote the solutions of © and m, respectively. An ensemble S(0) is called uniformly ensemble 
controllable, if for any x* € C(P,M”') and any e > 0 there exists a T > 0 and an input function 
u € L^([0, T], R”^) (or, in discrete-time, a finite input sequence uq, ...,ut-i) such that 

sup \\(p{T, 9,u) - x*{9)\\ < e. (3) 

eeP 

They key point to note is that the input function (or input sequence) is assumed to be independent 
of the parameter values 9. Thus the open loop input function u achieves ([3]) universally for all 9. For 
the subsequent analysis the following necessary conditions are important. For the proof we refer to [9l 
Lemma 1]. Let (t{A) C C denote the spectrum of an n x n-matrix A, i.e., the set of real and complex 
eigenvalues of A. 
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Proposition 1. Suppose that the ensemble m (or is uniform ensemble eontrollable. Then 
(El) For every 6 £ P the linear system {A{6),B{6)) is reachable. 

(E2) For each number s >m + l of distinct parameters 6i,... ,6s £ P, the spectra of A{9) satisfy 

a{A{ei)) n • • • n (T{A{es)) = 0. 

The preceding two conditions are useful in order to rule out families that are not ensemble con¬ 
trollable. In particular, condition (E2) implies that A{9) cannot have a ^-independent eigenvalue. We 
next show that necessary and sufficient conditions for uniform ensemble controllability can be stated 
in terms of a polynomial approximation property. 

For discrete-time ensembles with a scalar input sequence uq, the solution at T > 0 is 

ip{T, 9, u) = {uT-iI + ut-2A{9) -\ -h uoA{9)'^~^) b{9). 

This implies that a family of single-input, discrete-time systems {A{9),b{9)) is uniformly ensemble 
controllable if and only if for all e > 0 and all x* £ C'(P,M”') there is a real scalar polynomial p € M[ 2 ;] 
such that 

sup \\p{A{9)) b{9) — x*(0)|| < e. 

6»gp 

For a multivariable discrete-time system, the ensemble control condition ([3]) is similarly seen as being 
equivalent to 

m 

sup II ^pj{A{9)) bj{9) - x*{9)\\ < e. 

More generally, we obtain a polynomial characterization of ensemble controllability for both 
continuous-time and discrete-time systems. 

Remark 1. Let P be compact. Assume that: 

(Kl) A{9) has simple spectra for all 9. 

(K2) For all 9^9' the spectra of A{9)) and A{9') are disjoint. 

Then, the connected components of 


K=\J a{A{9)) C C 
6»eP 

are simply connected. Note that the converse is false in general. 

The next result shows that uniform ensemble controllability for continuous-time and discrete-time 
systems © and ([2]) are equivalent. 

Theorem 1. Let A £ and B £ Let bi{9),... ,bm{9) denote the columns of 

B{9). The following assertions are equivalent. 

(a) The continuous-time ensemble m is uniformly ensemble controllable. 

(h) The discrete-time ensemble ^ is uniformly ensemble controllable. 
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(c) For each e > 0 and each x* E there exist real scalar polynomials E M[z] 

such that 


m 

sup II '^pj{A{e)) bj{6) - x*{e)\\ < e. (4) 

,=i 


(d) The set 


span{A{9)^bj{9) \ k E No, j = 1,..., m} 

of continuous functions in 9 is dense in C{P,M^) with respect to the sup-norm. 

Proof. The preceding remarks prove the equivalence of (a) and (c). The equivalence of (c) and (d) is 
obvious and the equivalence of (b) and (d) follows from [171 Theorem 3.1.1], □ 


3 Refined characterizations of uniform ensemble controllability 


The preceding results show that the problem of uniform ensemble controllability is intimately con¬ 
nected to the classical problem of approximating a continuous function by a polynomial. In the sequel, 
unless stated otherwise, all results hold for uniform ensemble controllability of both continuous-time 
and discrete-time systems. In this section we will derive more explicit necessary and sufficient condi¬ 
tions. In particular, we will show that uniform controllability can only hold for systems depending on 
at most two parameters. 


We begin with a technical result. Given an array of linear parameter dependent systems {Aij , Bij ) E 
C'(P, X C'(P,M"'*^™'-’') with I < i < j < N. Let n = uud m = Define the 

associated upper triangular ensemble of systems by 


A = 

Mil • 

Aim ^ 

E G(P,R’^'''^), B = 

Mil • 

Bin\ 


1 0 

Ann / 


1 0 

Bnn ) 


(5) 


Proposition 2. The upper triangular family of systems {A{9),B{9)) is uniform ensemble controllable 
if and only if the families {Aii{9), Bii(9)) are uniform ensemble controllable for i = 1,..., N. 

Proof. The necessity part is obvious. Thus assume that {An, Bn) are uniform ensemble controllable. 
By Theorem d] it is sufficient to consider continuous-time systems. For simplicity we focus on = 2, 
i.e., on 


xi{t) = An{9)xi{t) Ai2{9)x2{t) + Bii{9)ui{t) + Bi2{9)u2{t) 

X2{t) = A22{9)x2{t) + B22{9)u2{t). 

The general case is treated, proceeding by induction. Let Xi{t) denote the solution of the i-th com¬ 
ponent. Given x* = col(a:]] • • • E C'(P,M”) and let e > 0. By uniform ensemble controllability of 
(^ 22 ,^ 22 ) there exists an input U 2 E L^([ 0 ,T],such that 

sup ||x 2 (r) - X2{9)\\ < S. 

0eP 

Let u = col(rii un) E L^([0,T],R™-) with Ui E L^([0,T],R™*). Applying the variations of 
constant formula we have ^ 

X2{t) = [ B22U2 {s) ds 

Jo 
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and thus 


xi{T) = zi{T)+ [ 

Jo 


where 


zi{T) = j ^^12 J dr + 512 ^ 2 ( 5 )^ ds. 

By uniform ensemble controllability of (^ii,i?ii) there exists an input ui G L^([0,T],R™'i) with 


sup II / ds — x^(0) + zi(T)|| < e. 

6»eP Jo 


But this implies 


and we are done. 


sup ||xi(r) - Xi(6»)|| < e 
eeP 


□ 


Using an appropriate similarity transformation, every system can be transformed into Hermite 
canonical form [ni DTO], which has the upper triangular form Q. Given a matrix pair {A,B) G 
•gnxn j^nxm^ where bi is the rth column of B. Select from left to right in the permuted Kalman 
matrix 

(61 Abi ■■■A^-% ••• bm Abm--- A^-^bm), 

the first linear independent columns. Then one obtains a list of basis vectors of the reachability 
subspace as 


bi,...,A^^-^bi,...,bm,...,A^^-^hm. 

The integers Ki, are called the Hermite indices, where Ki := 0 if the column bi has not been 
selected. One has Ki + • • • + Km = n if and only if {A,B) is reachable. The next result has been 
proven in [9] for the special case that P is a single compact interval. 

Theorem 2. Let P C M 6e a finite union of disjoint compact intervals. The ensemble of linear systems 
T, = {(A{9),B{6)) I 9 G P} is uniformly ensemble controllable if the following conditions are satisfied: 

(i) {A{9), B{9)) is reachable for all 9 £ P. 

(a) The input Hermite indices Ki{9 ),..., Km{9) of {A{9), B{9)) are independent of 9 £ P. 

(Hi) For any pair of distinct parameters 9,9' £ P, 9 9', the spectra of A{9) and A{9') are disjoint: 

a{A{9))na{A{9')) = ^. 

(iv) For each 9 £ P, the eigenvalues of A{9) have algebraic multiplicity one. 

Proof. We show the claim for P = Pi U P2, where Pi,P2 are disjoint compact intervals in M. Define 
Ki := UeeP '^iA{9)). The union of N disjoint compact intervals can be concluded by induction. Let 
X* £ C'(P,M"') and e > 0 be fixed. According to the proof of Theorem 1 in [9] it is sufficient to 
prove the assertion for discrete-time single-input ensembles. Then, by [9l Theorem 1] the ensembles 
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{{A{6),B{6)) \ 6 G Pi} and {{A{6),B{6)) \ 9 G P 2 } are uniformly ensemble controllable. Then, as is 
shown above there are polynomials pi G M[ 2 ;] and p 2 G M[z] such that 

sup ||pi(^( 6 ')) 6 ( 0 ) - x*( 6 ')|| < f and sup ||p 2 (^( 0 )) 6 ( 0 ) - x*( 0 )|| < 

6 »ePi 6 IGP 2 

Then, by the Stone-Weierstrass Theorem [ 6 l Theorem 6.6.3] there is a polynomial q G M[ 2 ;] satisfying 

sup \\pi{z) - g(z)|| < f and sup \\p 2 {z) - ^( 2)11 < |. 
zeifi z^K2 

Then, for any 0 G P, w.l.o.g. 0 G Pi, we have 

MA{9)) b{9) - x*(0)|| < ||g(^(0)) 6(0) - pM{0)) K^)ll + \\pi{A{9)) b{9) - x*(0)|| < e. 

This shows the assertion. □ 

In the continuous-time case, conditions (i)-(iv) imply that uniform ensemble controllability of S 
can be achieved in arbitrary time T > 0. As pointed out earlier, uniform ensemble controllability is 
related to approximation theory. The proof of Theorem 1 in is based on Mergelyan’s Theorem from 
complex approximation. 

Remark 2. Condition (iv) in Theorem\^ can be replaced by the following two conditions that are 
sometimes easier to cheek: 

1. A{9) is diagonalizahle by a similarity transformation with uniformly bounded eondition number. 

2. C\K is connected for K = IJeeP ^{A{9)). 

By Theorem [T] and Theorem 2.3 in [15j one obtaines for single-input ensembles the following 
sufficient condition. 

Theorem 3. Let P C K 6 e a compact interval and suppose the ensemble S = {(A(0),6(0)) | 0 G 
P} satisfies (i) and (in) in Theorem\M If there are ai,...,a„_i G M sueh that for all 9 G V the 
characteristic polynomials are of the form XA(e)(^) = z^ — an-iz"'~^ — ... — aiz — ao( 0 ), then S is 
uniformly ensemble controllable. 

As we next show, there do not exist ensembles of linear systems with more than three parameters 
that are uniformly ensemble controllable. We need the following lemma from |in| . 

Lemma 1. Let P C with int P = P. Let 11: —)■ 0 1 —). (A(0),R(0)) be real analytic. 

Then the Hermite indiees of the ensemble T, = {{A{9), B{9)) \ 9 G P} are generically constant, i.e. 
the Hermite indices are constant on an open and dense subset of P. 

Proof. Let Kn,m denote the set of all K = (Ki,..., Km) G N™ with ATi Km < n. For 

K G Kn,m let Her(Ar) := {{A,B) G x | Hermite indices are K} denote the subset in 

R"" of all systems with Hermite indices K. The set Her(Ar) is a constructible algebraic subset of 
R"" which induces a disjoint partition ^ Her(iF) = R” Consider the real analytic 

map H: R*^ —)■ 0 1 -), (A(0),i?(0)). Then, the preimage n“^(Her(Ar)) is a constructible 

analytic subset of R*^ (or the empty set) and there exists a AT* G Kn,m such that dimH”^ (Her(iF*)) 
has an interior point. But any constructible analytic subset 5 C R'^ of dimension d contains an 
open and dense subset of R'^. Thus the set of interior points of n“^( Her(iF*)) is an open and dense 
subset of R'^. In particular, the intersection H”^(Her(A'*)) n P contains an open and dense subset 
of P. It remains to show that iF* is unique. Suppose there exist two different Ki,K 2 G Kn^m such 
that their preimage H”^ (Her(Arj)), i = 1,2 contains open and dense subsets. Then the intersection 
n“^(Her(Ai)) n H”^ ( Her(Ar 2 )) 0, implying Ki = K 2 . This shows the assertion. □ 
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Theorem 4. Let P C with int P = P. Let 11: —>■ Q i->. (^A{9),B{9)) be real analytic. 

Suppose the ensemble S(0) = {[A{9),B{9)) \ 9 € P} is uniformly ensemble controllable. 

(a) Then dim P < 2. 

(b) Let Xi{9),..., Xn{9) denote the eigenvalues of A{9). Assume that at least one branch {Xk{9) \ 9 £ 
P} of the eigenvalues is contained in a one-dimensional real subspace S ofC. Then dim P = 1. 

(c) If there is at least one eigenvalue such that S := {Afc(0) | 9 G P} C M then dim P = 1. 

Proof. As the Hermite indices are generically constant, i.e. independent of the parameter, we may 
assume that {A{9),B{9)) satisfies the assumptions (i) and (ii) of Theorem [2j As the Hermite indices 
of the family {A{9),B[9)) are independent of parameter 0 G P, there exists a continuous family 
of invertible coordinate transformations S{9) such that {S{9)A{9)S{9)~^, S{9)B{9)) is in Hermite 
canonical form. Thus, w.l.o.g. we can assume that {A{9),B{9)) is in Hermite canonical form 


/An{9) • 

4im(0)\ 

T' ■ 

0 \ 

V 0 


Vo 

bm j 


where the m single-input subsystems {Ak}^{9),b}f) G xM”'' are reachable and in control canonical 

form. Note that b^ denotes the first standard basis vector and thus is independent of 9. Partition the 
desired state vector as x*{9) = col(x^(0) ••• x^(0)) with x^(0) G M"''. By Proposition [2] the single¬ 
input systems defined by {Akjf{9)^bjf) are uniform ensemble controllable. Applying (E2) for the single 
input case we conclude that the spectra of Akk{9) and Akk{9') are disjoint for all distinct parameters 
9^9' in P. Let A: P —>■ C be a locally defined branch of the eigenvalues of A[9). Let S denote 
be real subspace of C of smallest dimension such that A(P) C S. Since the eigenvalues of a matrix 
depends continuously on the parameters we see that locally A is continuous. By injectivity of A we 
get dim® P < dimjR S. This completes the proof. □ 


4 Special classes of ensembles 


Consider an ensemble of harmonic oscillators defined by 


9 A 




where the parameter is contained in the union P of hnitely many compact intervals with 0 G P. Note 
that for 0 = 0 the matrix A(0) has a double eigenvalue 0. The Hermite indices are 


Ki{e) 


2 if 0 / 0 
1 if 0 = 0 


and K2{9) 


0 if 0 / 0 
1 if 0 = 0, 


i.e. they are not constant. Thus the ensemble does not satisfy the condition (ii) in Theorem O 
Nevertheless we will show that this family of systems is uniformly ensemble controllable. The next 
result applies Theorem [2] to prove an extension of |14l Theorem 1]. 

For fixed matrices A G and B G we consider ensembles of the form 


S = {(0A,H)|0gP}. 

For this special class of system families the characterization of uniform ensemble controllability depends 
on whether the parameter set contains zero or not. 
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Theorem 5. Let P C M 6e the union of compact intervals. Assume that 0 G P. The family S = 
{{9A,B) I 9 G P} is uniformly ensemble controllable if and only if rank A = n and rankB = n. 

Proof. We focus on the continuous-time case; the discrete-time case goes mutatis mutandis. Suppose 
that S = {{9A,B) \ 9 G P} is uniformly ensemble controllable. Then, since 0 G P, the necessary 
condition (El) implies ranki? = n. In particular, we have m>n. To show the second claim, suppose 
that rank A < n. Then zero is an eigenvalue of A and for distinct parameter values {9i, ...,9n+i} G P 
we have 

0 G cr(6li^) n • • • n a(9n+iA) 
contradicting the necessary condition (E2). 

Conversely, assume that ranked = n and ranki? = n. Without loss of generality we can assume 
that B = In. The reachability condition (El) is implied by the rank condition on B. Let A denote 
the Jordan canonical form. It is sufficient to consider the ensemble 


dt 


x{t, 9) = 9Ax{t, 9) + I u{t). 


( 6 ) 


Using Proposition [2] it remains to prove the assertion for one Jordan block. Eor simplicity we focus 
on the case of a two dimensional Jordan block; the higher dimensional case follows by an induction 
argument. Let 


1 0 
0 1 


u{t). 


The solution to ([7]) is given by 

cT /^eX{T-s) _|_ 0 (7" _ 5 ) g 6 »A(T-s) .^^(s) 


(f{T,9,u) = 


^ex (T-s) 


U2is) 


ds. 


Given z* = col(zj[' € ^(P,!^^) and e > 0. By applying Theorem[2]to the ensemble 

^Z2{t,9) = 9Xz2{t,9) + u{t), 

there is an input function U 2 ■ [0, T] —>■ R so that 1^2 (^) “ T 2 {T, 9, U 2 )\ < s for all 0 G P. Let 

w*{9):=zl{9)- [ 0(r-s)e®^('^-")u2(s)ds G C(P,R). 

Jo 

Eollowing the same reasoning there is an input ui: [0, T] —)■ R so that 

rT 


W 


Consequently, we have 


and we are done. 


(9)- [ e^^(^-*)ui(s)ds| < e. 
Jo 


sup \\z*{9) — ip{T, 9, u)|| < e 
eeP 


(7) 


□ 


Theorem [5] dealt with the situation m>n. In the subsequent result we do not make this assump¬ 
tion. We use the notation AP := {A 0 | 0 G P}. 

Theorem 6. Let P be the union of compact real intervals with 0 0 P. 

(a) If the family T, = {{9A, B) | 0 G P} is uniformly ensemble controllable then (A, B) is controllable 
and A is invertible. 


(b) Let {A,B) be controllable and let A be invertible and diagonalizable such that AfcPnA;P = 0 for 
all k ^ I £ {1, Then the family S = {{6A,B) \ 6 £ P} is uniformly ensemble controllable. 

Proof, (a) Let the family S = {{OA, B)\6 £ P} be uniformly ensemble controllable. Then, by (El) the 
pair {6A,B) is controllable for every 0 € P. Using the Kalman matrix we have {A,B) is controllable. 
To show the second claim, suppose that rank A < n. Then zero is an eigenvalue of A and for distinct 
parameter values {0i, ...,0m+i} G P we have 

0 G cr(0i^) n • • • n a{9m+iA) 

contradicting the necessary condition (E2). 

(b) To show the claim, we verify the sufficient conditions of Theorem [2j The reachability of the 
pair {A,B) and the fact that 0 7^ P implies that {9A,B) is reachable for every 0 G P. Due to the 
fact that 0 0 P the Hermite indices of (9A, B) are independent of 0. Moreover, as A^P fl A;P = 0 for 
all A; 7^ / G {1, ...,r} we have a[9A) D a[9'A^ = 0 for all 0 7^ 0' G P. The assertion then follows from 
Remark O □ 

5 Ensembles of networks of linear dynamical systems 

We consider ensembles of networks of N identical single-input-single-output systems S = {A,b,c), 
whose dynamics are given by 


Xi{t) = Axift) + bvi{t) 

yi{t) = cxi{t). 

Here A £ b £ M"' and c G We assume that (H, b, c) is controllable and observable. The N 

identical systems S are coupled via directed links and the interconnection structure is described by a 
directed graph Q = (V, T). Here V = {1,... ,N} denotes the set of vertices, i.e. the N SISO systems, 
and £ describes the set of edges, i.e. the couplings. We examine the situation where the coupling 
strength is uncertain and is assume to vary over a compact interval P C M+ = (0, 00). The weighted 
graph adjacency matrix is given by 


m 


kij{9) ii{i,j)££ 

0 else. 


where kij : P —)■ (0, 00) are known continuous functions, {i,j) £ £. Thus, /Cp = {/C(0) | 0 G P} denotes 
a family of adjacency matrices. We assume that there is an external input u which is broadcasted 
to the systems S via the input-to-state interconnection vector B £ . The setting is illustrated in 

EigureO 

Thus, S = {(S, /C(0), H), 0 G P} defines an ensemble of networks. Using the state vector x = 
col(xi • • • xn) £ the network is described by the following control system 


■§ix{t, 9) = [l 0 A + /C(0) (g) bc)x{t, 9) + (^B 0 fe)u(t). (8) 

Let 1 = (1 • • • 1)"'" G and x^ £ M” denote the initial state of the ith node system and xq = 
col(X]’ • • • x^) £ The solution to ([8]) at time T > 0 starting in xg under the interconnection 

/C(0) and the input u is denoted by (p(T, /C(9), xg, u). 

We emphasize that the input u is broadcasted to the systems S within the network according to 
the input-to-state vector B. This phenomenon may be interpreted as that u serves as an universal 


9 


u 



Figure 1: Block diagram of an ensemble of homogenous networks with parameter-dependent couplings. 


input for a whole ensemble of networks that steers the initial state xq to a desired terminal state 
X* € in finite time T uniformly for all interconnection matrices JC{6), 0 G P. 

The ensemble {{T,,IC{9),B) \ 6 G P} is called robustly synchronizable to x* G C'(P,M"') if for every 
e > 0 there is a T > 0 and an inpnt-function u G L^([0,T],M) such that 

sup \\ip{T,K,{6),XQ,u) - (1 (8) x*(0))|| < e. 
e»GP 

The reason for this terminology is that the desired states to which we want to control have identical 
components, i.e., 

V(0)\ 


x*{e) = 


G 


T>nN 


yX*{9)J 


Thus the single input function achieves synchronization in finite time T, starting from non- 
synchronuous initial condition xq = col(xi(0), ■ ■ ■ ,X7v(0)) G We emphasize that the inputs 

may depend both on the initial condition xq and on the family of terminal state x*{9). 

As an illustration we consider N identical single-input-single-outpnt (SISO) harmonic oscillators 


A := 


0 -1 

1 0 


b := 


c :=[0 1 ], 


(9) 


The oscillators are coupled in a circular manner. The network topology is described by a directed 
graph Q with N nodes. The weighted adjacency matrix /C(9') is given by the circulant matrix 


^nng(9) : = 


0 9 


9 


0 GP := [ 9 -, 9 +]. 


( 10 ) 


The network is depicted in Figure [2j 

Without loss of generality, assume that the harmonic oscillators are numbered be such that the 
external input is applied to the first oscillator. Thus, the input-to-state interconnection vector is 
B = ei = (1,0, ...,0)"''. The dynamics of the overall network of systems is of the form 


^x{t, 9) = (/ (?) A -b JCring{9) (g) bc)x{t) + (ei (8> b)u{t) 
x{0,9) = xo G 


( 11 ) 
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Figure 2: Ensemble of rings of 5 identical harmonic oscillators. 

Let X* E denote the desired terminal states of the harmonic oscillators. 

Proposition 3. Let P denote a compact interval in (0,1) U (1, oo). The circular network ensemble of 
harmonic oscillators m is robustly synchronizable to 1 (8) x*. 


2 — 

Proof. Let oj := e N denote the primitive iV-th root of nnity. The adjacency matrix is a circnlant 
matrix with spectrum 


a{K{e)) = |/ = 0,...,iV-l 

The family {IC{0) \ 0 E P} of circulant matrices is simultaneously diagonalizable using the unitary 
Vandermonde matrix 


S = 


1 1 

1 OJ 


OJ 


1 

N-1 


1 UJ 


N-1 


OJ 


(Af-1)2 


( 12 ) 


Applying the change of coordinates S ® I to the network dynamics yields the state space eqnivalent 
system 

^z{t, 6) = diag {A + Xi{9)bc) x{t) + (1 (g) 6) u{t) 
z{0,0) = {S iS> I)xq. 


(13) 


Since S'! = Nei, the corresponding desired terminal states are Nei ® x*{9) E C'(P,M^'^). We show 
that Theorem [2] can be applied to (|13p . A simple calculation shows that 


N 


N 


y (T (A + \i{9)bc) = IJ J tc 


1=1 


1=1 


- {9e'^^^N - 1 ) = 0 


(14) 


Thus 


a{A + Xi{9)bc) n a{A + Xk{9)bc) = 0 

holds for all 0 E P and / 7 ^ A: E {1, Moreover, since {A + Xi{9)bc, b) is controllable for all 0 E P 

and I E {1,..., n}, we conclude that the parallel connected system (fT 3 l) is controllable. Furthermore, by 
inspection from (I14p . conditions (iii) and (iv) in Theorem[2]are satisfied. This completes the proof. □ 
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Proposition corrects an error in m, where the corresponding result was claimed for undirected 
circular graphs, with the symmetric circulant adjacency matrix 

'o e e 

^ . (15) 

e 

9 9 0 

This matrix has real eigenvalues 9 cos (27r^) , / = 0, N—1. Thus, for = 4, the necessary condition 
(E2) is not satisfied as cos(f) = cos(^) = 0. Therefore, the ring of oscillators with symmetric coupling 
matrix (1151) is not robust synchronizable. 

Next, we discuss a scenario where the node systems depend on a single parameter and are inter¬ 
connected by a fixed graph adjacency matrix JC. That is, let S(0) = {A{9),b{9),c{9)) be an ensemble 
of SISO systems, where A G C'(P,M”^”), h G C(P,M”) and c G C'(P,M^^’^). We assume that the 
system {A{9), b{9), c{9)) is controllable and observable for every 0 G P. Let Q = (V, T) be a directed 
graph with a fixed weighted adjacency matrix fC G 



Figure 3: Block diagram of the ensemble of homogenous networks with parameter-dependent node 
systems. 


Apply a single, parameter-independent external input u to the systems S(0), using an input-to- 
state interconnection vector B G Let x = col(xi • • • xn) G denote the sate vector of the 
network. The ensemble of networks we are interested in is described by the control system 

|x(L 0) = (/ 0 A{9) + /C 0 b{9)ci9))x{t, 9) + {B0 6(0))u(t). (16) 

We note that (fTHl) is a special case of the more general setting described in [8]. In particular, the 
controllability properties of the network are easily established. In fact, (1161) is controllable if and only 
if {IC,B) and {A{9), b{9), c{9)) is controllable for all 9] cf. [8]. In order to simplify the analysis, we 
assume that the adjacency matrix 1C has distinct eigenvalues Ai,..., Atv, which implies that there is an 
invertible matrix S with S1CS~^ = V = diag(Ai,..., Xn)- Using the change of coordinates z = (S'0/)x, 
system dUD is state-space equivalent to 


f z(L 0) = (/ 0 A{9) + P 0 b{9)c{9))z{t, 9) + {SB 0 b{9))u{t). 
Note that A^ ^ A; for all A: 7 ^ / G {1,...., A"}. System (fT7)l is equivalent to 
(A{9) + Xib{9)c{9) 


( 17 ) 


§-Mt,9) = 


z{t, 9) -|- {SB 0 b)u{t). 


V 


A{9) + XNb{9)c{9). 
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Thus the spectrum of I ^ A{6) + /C 0 b{9)c{6) is 

a {A{e) + Ai b{e)c{e)) u • • • u u {A{e) + Xn b{e)c{e)). 

Since {A{9),b{9), c[9)) is controllable and observable we obtain 

det [zl - {A{9) - Xkh{9)c{9))) = qg{z) + XkPe{z) 

with coprime polynomials peiz), qe{z). This shows, for fixed 0 G P and Xk ^ Xi, that there is no z G C 
such that 


qe{z) + XkPeiz) = qg{z) + Xipe{z). 

Therefore, one can apply Theorem [21 which proves the next result. 

Proposition 4. Let P C M 6e compact. Suppose the system S(0) = {A{9), b{9), c{9)) is controllable 
and observable for every 0 G P. Assume that the pair {IC,B) G is controllable and 1C is 

diagonalizable. Then the ensemble of networks (HID is robustly synchronizahle if the eigenvalues of 
A{9) + A b{9)c{9) have algebraic multiplicity one for each 0 G P and X G ct(/C) and for each 9 9' 

and Xi ^ Xk € (^{IC) it holds 

a {A{9) + Xi b{9)c{9)) n u {A{9') + Xk b{9')c{9')) = 0. 


6 Conclusions 

An approximation theoretical approach to controlling state ensembles of linear systems has been 
proposed. We concluded that the ensemble control problem for continuous-time and discrete-time 
systems is equivalent. For analytic parameter-dependent linear systems it is shown that uniform 
ensemble controllability is not possible for more than three parameters. A complete characterization 
of uniform ensemble controllability is presented for a special class of ensembles. An application to 
robust synchronization using broadcasted open-loop controls is given. Our approach is based on 
information of the spectrum of the system matrices (see assumption (iv) in Theorem [2]). An open 
problem is to find relaxed necessary and sufficient conditions using the concept of pseudospectra. 
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